Two-dimensional topological solitons in rectangular magnetic dots. 



Konstantin L. Metlov 
Institute of Physics ASCR, Na Slovance 2, Prague 8, CZ-18221 
(February 1, 2008) 



O 
O 

(N 



I 

-a 
c 

o 
o 



> 

o 
in 
o 



i 

C 

O 

o 



A general approach allowing to find the analytical ex- 
pressions for equilibrium magnetic structures in small and 
flat magnetic nano-sized cylinders of arbitrary shape made of 
soft magnetic material is presented. The resulting magneti- 
zation distributions are two-dimensional topological solitons 
and have a non-zero topological charge. The approach is il- 
lustrated here on an example of a thin rectangular particle. 
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Introduction. The small magnetic cylindrical ele- 
ments (dots) made of a soft magnetic material became 
the focus of attention due to their potential applications 
in the Magnetic Random Access Memory (MRAM) de- 
vices. These elements can be made of different shapes 
and sizes (usually of the order of 10 — lOOnm ) and 
found to display a variety of magnetic structures (equi- 
librium distributions of magnetization vector field). The 
structures are solutions of non-linear and non-local (inte- 
gral partial differential) equations of micromagnetics and 
were a recent subject of extensive study by finite element 
computer simulations (see e.g. and references therein). 

Recently, a new method allowing to treat the problem 
of finding the magnetic structures of thin dots analyti- 
cally became available |^|| . It is based on the existence 
of the well defined hierarchy of energies in small flat mag- 
netic cylinders made of soft magnetic material || . The 
resulting analytical magnetization distributions minimize 
the exchange energy exactly, have no magnetic charges 
on the cylinder faces (except in a few points of topologi- 
cal singularities) and no magnetic charges on the cylinder 
sides. 

Briefly, the general recipe of Ref. || can be formu- 
lated as follows. In a thin ferromagnetic cylinder (with 
thickness L < Le, where Le — y/C/Ms is the ex- 
change length, C is the exchange constant, Ms is the 
saturation magnetization of the material) the dependence 
of the magnetization distribution rh{r) = M(r)/Ms on 
the coordinate Z perpendicular to the cylinder face can 
be neglected (so that m(X,Y), where X and Y are 
Cartesian coordinates on the planes parallel to the cylin- 
der face). Also, because the magnetization vector has 
a fixed length |m| = 1 it has only two independent 
components. Both these facts make it convenient to 
parametrize the magnetization vector field by a com- 
plex function w(z) of the complex variable z = X + lY, 
i = \J — 1, so that mz = (1 — w(z)w(z)) / (1 4- w{z)w{z)), 
mx + imy = 2w(z)/(l + w(z)w(z)), where a line over 



symbol means the complex conjugation. The require- 
ment of minimum face magnetic charges makes the func- 
tion w(z) non-analytical (in the sense of the Cauchy- 
Riemann conditions), but it is possible to express it 0] 
through another, analytical, function f{z). For a particle 
with corners this can be done [H] as 



/(*) I/C0I < 1 

w ( z ) = { f(z)/yff(z)J(zj K\f(z)\<d 
f(z)/d > d, 



(1) 



where d is an arbitrary constant defining the size of re- 
gions of the out-of-plane magnetization vector at topo- 
logical singularities at cylinder corners. 

Let us denote the conformal transformation of the in- 
terior of the unit circle \t\ < 1 to z € T> (T> is the set 
of points of the X — Y plane belonging to the cylinder 
face) as z = T(t), according to Riemann theorem such 
transformation exists for any simply connected region T>. 
Then, the analytical function f(z) corresponding to the 
magnetization vector field with no magnetic charges on 
the cylinder side is given || by 

f(z) = (itc + A- At 2 )T'(t), t^T-\z), (2) 

where c and A are arbitrary real and complex constants 
respectively, T^ 1 (z) is the inverse function of T(t). The 
function T~ 1 (z) exists and is unique because the trans- 
formation is conformal. The expression (|^) contains es- 
scntilly all the solutions minimizing the exchange energy 
and having no magnetic charge on the cylinder sides. 
The magnetostatic energy of the resulting solutions is 
also very small compared to the exchange energy (but 
this depends on the particle size) and it can be expected 
(for sufficiently small particles) not to introduce any new 
solution types but just to modify slightly. 

a. The magnetic structures in a rectangular particle. 
The shape of the particle enters the resulting distribution 
through the conformal transformation T(t). The confor- 
mal transformation of the unit circle to the rectangle is 



T(t)=d( 



o VT + t 4 - 2t 2 cos 26 



dt + Co), 



(3) 



which is a particular case of Schwartz-Christoffel integral. 
This integral can be expressed analytically through the 
elliptic integral of the first kind F(<fi, m). The parameter 
< 5 < 7r/4 controls the aspect ratio of the rectangle, 
<5 = 7r/4 corresponds to square, S ~ 0.172426 to the 2 : 1 
rectangle. The constants C\ and Co were chosen in such 
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a way that the height of the rectangle is 1 and its upper 
left corner is situated at the point with X, Y = 0. For the 
2 : 1 rectangle shown in the pictures C\ ~ 0.63189 and 
Co — —1.849 + i0.787. This transformation is visualized 
in Fig. |. 

The magnetic structures in the particle are thus 



f(z) = (itc + A- At 2 ) 



VI + t 4 - 2t 2 cos 25 



t->T~\z), 
(4) 



where C\ is absorbed into c and A 

To find the values of constants c and A it is required 
to minimize the total energy of the particle (including 
magnetostatics). This will be done elsewhere, let us now 
just treat these constants as independent parameters and 
analyze the possible types of solutions qualitatively. 

As it was already shown || for the cases of circular 
and triangular cylinders (and it is applicable to any soli- 
ton of type (||)) there are two basic classes of solutions 
depending on whether \c\ < 2\A\ or not. The first class, 
corresponds to the case |c| > 2\A\ is a vortex, the sign 
of c controls the direction of the rotation of the magneti- 
zation vector in it. If the particle is square the centered 
vortex (Fig ||a) is stable. In elongated particles a state 
with displaced vortex (Fig||b, see also Figs 3,5 in Ref. 
can also be. The other class of solutions contained in (0) 
corresponds to the two skyrmions (hedgehogs) bound to 
the cylinder sides. If the particle is not of the circular 
shape positions of the hedgehogs on its boundary are not 
equivalent. Particularly, for particle with corners the to- 
tal energy is lower if the hedgehogs are located in the cor- 
ners (however states with hedgehogs at the middle of the 
edges can be metastable). For rectangle (because of the 
symmetry) there are two possibilities: either the hedge- 
hogs are located at the neighbouring corners (Fig||c), or 
at the opposite ones (Fig ||d) . Both these structures, usu- 
ally called "C" and "S" , already appeared in numerical 
micromagnetic simulations for rectangles. 

Because of the completeness of the solution (||) no 
other shapes of the equilibrium magnetization vector field 
can be expected in the rectangular particle if it is suf- 
ficiently small (but larger then the single domain size, 
where absence of magnetic charges is not required) and 
if it is made of sufficiently isotropic magnetic material, 
which are the basic assumptions of the presented theory. 

The expression ([|) also describes the metastable states 
when the c and A deviate from their equilibrium values. 
These states can be useful for considering the magneti- 
zation reversal, excitations and the stability of the equi- 
librium solutions (see e.g. 0). 

Summary. A simple analytical framework for find- 
ing the profiles of two-dimensional topological solitons in 
small flat ferromagnetic cylinders of the arbitrary shape 
was presented. As an example, the case of the rectangu- 
lar cylinder was considered. 



The resulting magnetization vector fields may find 
their use in magnetic microscopy (for fitting the signal 
from the probe by changing the parameters A and c), 
and for building the phase diagrams of small magnetic 
particles outlining the regions of the particle geometri- 
cal parameters where each of the presented equilibrium 
solutions correspond to the ground state of the system. 
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FIG. 1. Conformal mapping of the unit circle to the rect- 
angle, dots on both plots correspond to each other to show 
how the interior is transformed. 
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FIG. 2. The equilibrium configurations of the magnetiza- 
tion vector (its projection to the cylinder face is shown) in 
the rectangular particle (|l]), (jjj) with d = 6: a) vortex c = 2, 
A — 0, ; b) displaced vortex c = 2, A — 1.9i; c) "C" structure 
c = 2, A = 11.6; d) "S" structure c = 0, A = 6e'° ' 17242B . 
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